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ABSTRACT 

Hypercritical accretion flows onto stellar mass black holes (BHs) are com- 
monly considered as a promising model of central engines of gamma-ray bursts 
(GRBs). In this model a certain fraction of gravitational binding energy of accret- 
ing matter is deposited to the energy of relativistic jets via neutrino annihilation 
and/or magnetic fields. However, some recent studies have indicated that the 
energy deposition rate by neutrino annihilation is somewhat smaller than that 
needed to power a GRB. To overcome this difficulty, Ramirez- Ruiz & Socrates 
(2005) proposed that high energy neutrinos from hot corona above the accretion 
disk might enhance the efficiency of energy deposition. We elucidate the disk 
corona model in the context of hypercritical accretion flows. From the energy 
balance in the disk and the corona, we can calculate the disk and coronal temper- 
ature, Td and T c , and neutrino spectra, taking into account the neutrino cooling 
processes by neutrino-electron scatterings and neutrino pair productions. The 
calculated neutrino spectra consist of two peaks; one by the neutrino emission 
from the disk and the other by that from the corona. We find that the disk 
corona can enhance the efficiency of energy release but only by a factor of 1.5 
or so, unless the height of the corona is very small, H <C r. This is because the 
neutrino emission is very sensitive to the temperature of emitting region, and 
then the ratio T c /Td cannot be so large. 

Subject headings: accretion, accretion disks — black hole physics — gamma rays: 
bursts — neutrinos 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most explosive phenomena in the Universe. Recent 
development in the observational study of GRBs and their afterglows lead to the under- 
standing of their emission mechanisms (Piran 2005; Meszaros 2006 for reviews). The prompt 
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emissions of GRBs are thought to be emitted from optically thin plasmas with highly rela- 
tivistic velocity in order to avoid the compactness problem (Krolik & Pier 1991; Lithwick & 
Sari 2001). 

Some observational studies show jet breaks in their lightcurves, implying the relativistic 
outflow may be anisotropic, or jet (Rhoads 1999; Sari, Piran, & Halpern 1999). GRB prompt 
emission has an energy E 1 ~ 10 51 ergs, if we take into account the beaming of radiations 
(Frail et al. 2001; Bloom, Frail, & Kulkarni 2003), and the kinetic energy of relativistic jet 
producing GRB is about E tot ~ 10 52 ergs, if we assume the conversion efficiencies of kinetic 
energy into radiative energy as about 0.1 (Beloborodov 2000; Guetta, Spada, & Waxman 
2001). 

The engines powering such energetic outflows are still unknown. Some long GRBs 
show the associations with supernovae type Ic, implying that they originate from the core 
collapse of massive stars (Bloom et al. 1999; Woosley 1993). Short GRBs are, on the other 
hand, recently found to be associated with elliptical galaxies (Gal- Yam et al. 2005), so they 
may originate from neutron star (NS) mergers or black hole(BH)-NS mergers (Paczynski 
1991; Narayan, Paczynski, & Piran 1992). As the result of such energetic phenomena the 
gravitational energy released in those processes can drive the relativistic jets. A plausible 
model of the central engine of GRBs involves an accretion disk with extremely high accretion 
rate (M ~ 0.01 — 1OM ) around a stellar mass BH. Such a flow is called as a neutrino cooled 
accretion flow or a hypercritical accretion flow (Kohri & Mineshige 2002). Recent studies of 
hypercritical accretion flows have been developed by Narayan, Paczynski, & Piran (1992), 
Popham, Woosley, & Fryer (1999), Di Matteo, Perna, & Narayan (2002, hereafter DPN), 
Kohri, Narayan, & Piran (2005) and Kawanaka & Mineshige (2007). The radiative cooling 
is inefficient in such a hot, dense accretion disk with hypercritical accretion rate since the 
optical depth of disk is so high that photons cannot escape thorough the disk within the 
accretion timescale. The dominant cooling process of hypercritical accretion flows might be 
neutrino processes with weak interaction and the total energy of neutrino emission reaches 
some fraction of the rest mass energy of accreting matter (DPN). The released gravitational 
energy as neutrino emission can deposit its energy to the relativistic fireball via neutrino- 
antineutrino annihilation in the baryon poor region along the rotational axis. However some 
previous work indicated that the energy released via neutrino process might be insufficient 
for GRB. The available energy for deposition by neutrino annihilation is only a few tenths 
of a percent of total energy of emitted neutrinos (Ruffert et al. 1997). DPN also pointed 
out that energy deposition rate by neutrino process might be inefficient since the accretion 
flow with a high accretion rate becomes thick for neutrinos and generated neutrinos cannot 
escape from accreting matter (so called as neutrino trapping). 
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Ramirez-Ruiz & Socrates (2005, hereafter RS) indicated that non-thermal neutrinos 
generated in the hot corona (with temperature T c ) above the hypercritical accretion disk 
(with temperature T^) might enhance the energy deposition by neutrino annihilation since 
the deposition rate is proportional to the mean energy of neutrinos and T c can be greater 
than Td. They assumed the coronal depth for neutrinos to estimate the coronal temperature 
and concluded that energy deposition rate may increase by some factor if a sufficiently thin 
corona forms. However they did not calculate the neutrino spectra and mean energy of 
emitted neutrinos. Moreover the coronal region they assumed may be too small, and thus, 
the coronal temperature might have been overestimated. 

In this paper we elucidate the disk corona model in the context of hypercritical accretion 
scenario and show that the mean energy of emergent neutrinos from hypercritical accretion 
flows is enhanced by the existence of neutrino-thin corona. We take a reasonable value of 
the coronal thickness and calculate the emergent spectra of neutrinos from the corona. In 
§2 our model of disk corona is described with some adopted assumptions. We show the 
numerical method to calculate the neutrino spectra, and energy equations to determine the 
temperature of the corona and the disk are written in §3. We consider the absorbed neutrinos 
by the dense disk, which re-emits thermal neutrinos, to determine the coronal temperature 
self-consistently, which RS did not take into account. Then we present our results in §4 and 
discuss the possibility of the corona to enhance the energy deposition rate in §5. 



2. OUR MODEL 

We first outline our model in this section. Originally, disk corona model was devel- 
oped to explain the non-thermal spectral component of X-ray binaries and active galactic 
nuclei (AGNs) in analogy with the model for the solar corona (see e.g. Shibata, Tajima, & 
Matsumoto 1990). From the fitting to the observed spectra it has become clear that /, the 
fraction of gravitational energy released in diffuse corona, should be almost unity, i.e, almost 
all the liberated gravitational energy should be dissipated in the corona, since otherwise the 
large flux of the non-thermal spectral component, which is compatible to that of the thermal 
component, cannot be explained (Haardt & Maraschi 1991). The magnetic pressure can be 
amplified up to the equipartition value with the gas pressure by MHD processes, such as 
magnetorotational instability in the accretion disk (Balbus & Hawley 1998; Machida & Mat- 
sumoto 2003). Magnetic loops are formed and lifted up by the Parker instability, creating a 
corona filled with magnetic fields (Galeev, Rosner, & Vaiana 1979; Stella, L., & Rosner, R. 
1984). This corona may be heated by magnetic reconnection. In this way some fraction / of 
total gravitational energy can be released in the corona. Liu, Mineshige, & Shibata (2002) 
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have shown that / ~ 1 from the consideration of energy balance in the disk corona system. 

Hot corona might also form in hypercritical accretion flows since the accretion disk is 
highly turbulent (RS). Thus, following the ordinary disk corona model mentioned above, we 
construct the disk corona model in the context of hypercritical accretion flows for GRBs. The 
high energy neutrinos from the hot corona can enhance the energy deposition rate by neutrino 
annihilation, since the annihilation cross section depends on the energy of neutrinos. The 
energy deposition rate by neutrino annihilation at per unit volume is (Ruffert et al. 1997) 

d- e -*f = 2(Cv2 3 + Ca2) <v J dn v J dn P ((e u ) + (e p )) - cosfl) 2 , (i) 

where C\ = 1/2 + 2sin 2 0\y and Ca = 1/2 for electron type pairs, Cy = —1/2 + 2sin 2 6 l w 
and Ca = — 1/2 for heavy-lepton neutrino pairs, Gp is the Fermi coupling constant, (E u ) = 
J EI v {E)dE/ J I u {E)dE (or (E p ) = J EI p {E)dE/ J I p {E)dE) is the mean energy of neu- 
trinos (antineutrinos), l v (or I p ) is intensity of emitted neutrinos (antineutrinos), and 9 is 
the collision angle, and VL U (or Q p ) is the solid angle of neutrino (antineutrino) emission, 
which depend on the geometry of emitting region. We take the Weinberg angle, 9-w, as 
sin 2 #w = 0.23. For simplicity, we assume that neutrinos and antineutrinos are all electron 
type and have the same intensity. Thus the energy deposition rate is roughly proportional 
to the mean energy of neutrinos. If the hot corona emitting neutrinos form, the neutrino 
spectra, which are generally assumed as thermal, may be deformed and have higher mean 
energy than original value. 

The energy released in the accretion disk per unit surface can be calculated by using 
the standard disk theory and is written as (Kato, Fukue, & Mineshige 1998) 

87rr d \ V r / 

where M, M, and r in is the mass of black hole, the accretion rate, and the inner boundary 
radius of the accretion disks where the torque vanishes, respectively. We take M = 3M Q 
and r; n = 3r s ~ 3 x 10 6 cm where r s is the Schwarzschild radius. This dissipated energy 
is, if the neutrino cooling is efficient, released as neutrino emission. We consider the plane- 
parallel and homogeneous corona with vertical thickness H above the disk (Fig{T]). With 
conservative views the parameter value H may be comparable to the disk thickness, since 
the most unstable wavelength of Parker instability is nearly scale height in the disk (e.g. 
Matsumoto et al. 1988) and thus, the typical scale of the magnetic field in the corona would 
be comparable to the scale height. 

Neutrino emission from the inner part of the disk dominates over that from the outer 
part and mainly contributes to the heating of relativistic fireball above the disk. Then, 
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we evaluate the enhancement of mean energy of neutrinos at r = 4r s , where the energy 
dissipation rate reaches nearly a maximum value, and we assume that this enhancement is 
proportional to that of the total energy deposition rate by neutrino annihilation. We assume 
that the corona consists of pure relativistic electron-positron plasma since the corona may 
form above the surface of the disk where the density of baryon is much less than that inside 
of the disk (discussed later). We also assume that the electrons and positrons are completely 
thermalized, since the timescale of electromagnetic interactions is much shorter than that of 
weak interactions by many order of magnitude. 

As the coronal cooling processes we take into account neutrino reactions as 

v + e ± -> v + e ± (3) 
+ e + — >■ v + v. (4) 

The thermal neutrinos emitted from the disk are up-scattered by hot electrons and positrons 
in the corona. The scattered neutrinos have higher energy since the temperature of the 
corona (T c ) is higher than that of the disk (T d ) and, hence, emerged neutrino spectrum is 
deformed. The corona is cooled via scatterings of neutrinos by high energy electrons and 
positrons in the corona. Some part of scattered neutrinos are re-absorbed and heat the disk, 
which re-emits thermal neutrinos. Neutrino pairs produced by annihilation of high energy 
electrons and positrons in the corona also cool the corona and have high energy spectra. 
We neglected the scattering of neutrinos produced by pair process in the corona since the 
energy of neutrinos emitted in the corona is almost the same energy as that of hot electrons 
and since the optical depth of electron-neutrino scattering is much less than unity in the 
corona. We also neglected the Fermi blocking effect by background of electrons, positrons 
and neutrinos. 

On the basis of these assumptions, we can calculate T d and T c by solving energy equa- 
tions including neutrino processes, which are needed to calculate the neutrino spectra (see 
§3.3). 

3. CALCULATIONS OF NEUTRINO SPECTRA 

Hereafter we take the units as c—1, h — 1, and k—1. Calculating the neutrino spectrum 
emitted from the disk corona system is important when we evaluate the energy deposition 
rate by neutrino annihilation, since the deposition rate is proportional to the neutrino energy. 
In this section we present the method to calculate the coronal neutrino spectra and their 
mean energy. 
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3.1. Neutrino-electron Scattering 

We calculated the change of energy of neutrinos up-scattered by hot coronal electrons 
and positrons. The differential cross section of neutrino-electron scattering is ('t Hooft 1971; 
Kneller, McLaughlin, & Surman 2006) 

da ev G F 2 m e 2 [ t „ f ^ 2 2 K e 2 



6K. 2, m , ^ + c t y + { c v -c A m-^y 

-(C V 2 -C A 2 )^>, (5) 

where Cy = 1/2 + 2sin 2 6 l w and Ca = 1/2 for electron-electron neutrino scattering, Cy = 
1/2 + 2sin 2 w and C A = —1/2 for positron-electron neutrino scattering, K e is the kinetic 
energy of recoiled electron, and E u is the energy of an incident neutrino in the rest frame of 
electron. We do not deal here with the composition of neutrino flavors, since the detailed 
composition depends on the structure of the disk. We assume that the emitted neutrinos 
from the disk are all electron type since nucleon pair capture is a dominant neutrino process 
in the disk. 

The mean free path, A, of neutrinos in the corona for neutrino-electron scattering is 
then (Landau & Lifshits 1976) 



A {2nf J Ml-/i/?)/A (6) 

where /3, p, f e are the velocity, momentum, and the distribution function of electrons, re- 
spectively, and ji is cosine of the angle between the direction of neutrino and electron. We 
assume that the corona consists of pure electron positron plasma in thermal equilibrium by 
electromagnetic process so that the distribution function of pair should be Fermi-Dirac type, 
f e = [exp(y / p 2 + mf/T) + 1] _1 . In the relativistic limit, where T c ^> m , the total coronal 
depth for neutrino-electron scattering is 



H 

T " A 

°' 27 ( v 10 6 cmJ V^lOMevJ V^lOMevJ ' ^ 

where E v is the energy of incident neutrino in the laboratory frame. Hence, higher energy 
neutrinos have more chances to collide with coronal electrons than lower ones. 



H \ ( E„ \ ( T c V 



We calculate the neutrino spectra scattered by coronal electrons and positrons by using 
Monte Carlo method, following Pozdnyakov, Sobol, & Sunyaev (1977) and Liu, Mineshige, 
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k, Ohsuga (2003). We first set the weight of neutrino w = 1 for a given thermal neutrino 
with energy, E v , which has the Fermi-Dirac distribution. We calculate the probability for 
neutrinos passing through the corona, Pq = exp(— rj cos a), where a is the angle between 
neutrino direction and the z-axis which is perpendicular to the coronal plain. Then WqPq 
is the transmitted portion of neutrinos and the remaining w\ = wq(1 — Pq) is the portion 
of neutrinos scattered at least once. Let w n = u> n -i(l — P n -i) be the portion of neutrinos 
experiencing the n-th scattering. We continue the calculation until the weight, w n , becomes 
sufficiently small. Repeating the same procedures for sufficiently large number of neutrinos, 
we can calculate emergent spectra by collecting neutrinos going upward through the coronal 
surface at z — H, while downward neutrinos crossing the lower boundary are re-absorbed 
by the disk body, thereby heating the disk. 



3.2. Neutrino-antineutrino Pair Production 

We also calculate the neutrino pair emission by weak interaction in the pair plasma. 
The cross section for pair process by weak interaction is (Dicus 1972; Yakovlev et al. 2001) 

G F 2 m c 4 



V = 



12lT E\E 2 



+3(CV - CV) 1 + 2 



m e 2 



(8) 



where Cy = 1/2 + 2sin 2 6 l w and Ca = 1/2 for pair production of electron types, Cy = 
— 1/2 + 2sin 2 # w and Ca = —1/2 for pair of /i and r neutrinos, Si (or E 2 ) and Pi (or P 2 ) 
are the energy and four momentum of electron (positron), respectively, and v is the relative 
velocity of the pair. The neutrino pair emissivity is given as 

g e -e+-i/e = ^2~~p j d 3 pid 3 p 2 (E 1 + E 2 )a e - e +^ l/P vf 1 f 2 . (9) 

We assume that the distribution functions of electrons and positrons are Fermi-Dirac type 
with the same temperature, T c . In the relativistic limit where T c ^> m e , the total emissivity 
becomes 

q e - e+ ^ = 1.39 x 10 34 (y^f) ergs cm" 3 s" 1 (10) 

and the mean energy of neutrinos is 

(E) j ' d z v x d z y 2 E^a vv vf x f 2 
J d*p l d?p 2 E l cr vv vf l f 2 
~ 5.1 T c . (11) 
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A half of the emitted neutrinos pass through the corona upwardly and the remaining half 
are absorbed and heat the disk since emission of pair process is isotropic. The emissivity, 
equation (jTUJ) includes all types of neutrinos but we assume that the emitted neutrinos 
by pair-process are electron type only when considering heating of the disk by downward 
neutrinos and the energy deposition to the relativistic fireball by upward neutrinos. This 
is a reasonable assumption in doing simple estimation, since the ratios of the emissivity of 
each type of neutrinos is, from equation flU]), 0.70:0.15:0.15 for u e , v T . 

We also calculated the neutrino spectra from pair process by Monte Carlo method. The 
reaction rate of electron positron annihilation by weak interaction is given as 

4 

dR e - e +^ v9 = -j—^a e - e +^ vV v fid 3 p 1 f 2 d 3 p 2 . (12) 

Thus, for given energy of electrons and positrons the reaction weight, w, is proportional to 
a e - e +^ u pV . The differential cross section of neutrino pair process in the frame of the center 
of mass is (Misiaszek, Odrzywotek, & Kutschera 2006), 

da -~- + ^ oc \ M \ 2 
dVt 

= 8G F 2 [(C v - C A ) 2 P 1 ■ Q 1 P 2 ■ Q 2 

+ (C Y + C A ) 2 P 2 ■ Q 1 P 1 ■ Q 2 + m e 2 (C v 2 - C A 2 )Qx • Q 2 ] , (13) 

where Pi (or P 2 ) and Qi(or Q 2 ) are the four momentum of electron (positron) and neutrino 
(antineutrino), respectively, and \Ai\ 2 is the amplitude of pair process. We determine the 
emitted energy of neutrino pair with reaction weight w from the differential cross section 
formula and calculate the neutrino spectra from pair process. 



3.3. Energy Balances in the Corona and the Disk 

We calculate the temperature of the corona and that of the disk consistently by solving 
their energy balances simultaneously and by calculating the neutrino spectra (see Figf2]). 
We assume that the fraction / of the total energy is dissipated in the corona and remaining 
part, 1 — f, in the disk, neglecting the advection of energy. The energy balance in the disk 
is 

(1 - f)Q+ + Q ref + \q e - e+ ^ uD H = 7 -aT d \ (14) 

where we assumed that the corona is thin for neutrinos emitted by pair process so that 
we can neglect the absorptions by the inverse pair process. On the other hand, the energy 



- 9- 



balance in the corona is 

7 

fQ + + oCiTd 4 = Q esc + Q rc{ + q e - e +_„ p H, (15) 
o 

where Q esc is the cooling rate by the thermal neutrinos from the disk and upward neutrinos 
scattered by electrons and positrons, and Q re f is that by the downwardly scattered neutrinos, 
which can be evaluated by calculating the scattered neutrino spectra. 

If all of the energy is dissipated within the disk, i.e., when f=0, the effective temperature 
of the disk at r = 4r s is 

/ M \ ~ 1/2 / M \ 1/4 

r * = 4J W \Tm^) MeV ' (16) 

and the mean energy of thermal neutrinos from the disk is (E)q ~ 4.1T<i,o- With given 
parameters, / and H, we can solve these equations ffT4l) and (ITB1) iteratively for T& and T c , 
together with the calculated neutrino spectra. Then, we evaluate the enhancement of the 
mean energy of neutrinos from the corona, comparing with that of thermal one, A.IT^q. 



4. RESULTS 

4.1. Neutrino-electron Scattering 

First we consider the effect of scattering alone as a cooling mechanism of the corona 
in analogy with the coronae of standard accretion disks in X-ray binaries and AGNs. From 
equations (JTlj) and (}Po]) . and neglecting the pair process, we can calculate the temperature of 
the corona and the disk by calculating the change of neutrino energy Q esc and Q Ie { by Monte 
Carlo method iteratively. Figure [3] shows the coronal and disk temperatures for M = lM s _1 
as functions of /. A thinner corona (with smaller r) has a higher temperature for a given 
Td, since the cooling rate of the corona by scatterings is roughly written as 

fQ+ ~ 7 -aT d ' x r x |f , (17) 

o Id 

where we assume that the scattered neutrino has nearly the same energy of that of electrons, 
and that r is proportional to the thickness of the corona. 

We then show the neutrino spectra where / is close to unity (FigHJ left) . The neutrino 
spectra have two components: thermal neutrinos from the disk and scattered neutrinos in the 
hot corona. Nearly a half of the scattered neutrinos go downward and are absorbed by the 



-10- 



disk, which generate soft thermal neutrinos, and the remaining half pass upward producing 
a high energy spectrum. Once thermal neutrinos are scattered, they acquire energy nearly 
equal to that of the coronal electrons, since the energy of neutrino is too high (T^ 3> m e ) to 
be scattered elastically. The high energy part of emergent neutrino spectra mainly contribute 
to neutrino heating by neutrino annihilation. If H is on the same order of r, in contrast, the 
the coronal temperature is not significantly high compared with disk temperature. Thus, 
the effect of up-scatterings is small and the spectrum does not show two humps so clearly. 

We also calculate the mean energy of neutrinos from neutrino spectra and the ampli- 
fication of neutrino mean energy is shown in Figure H] (right). With conservative value of 
coronal thickness, H/r ~ 0.1 — 1, the amplification factor of mean energy is about a factor 
of two. 



4.2. Neutrino Pair Production 

In the next case we take into account the effect of cooling by neutrino pair produc- 
tion. Neglecting the scattering, the temperature of the corona and the disk can be solved 
analytically. The temperature of the corona and the disk at r = 4r s are 

-1/9 

MeV (18) 

1/4 

MeV, (19) 

and the mean energy of emergent neutrinos is (E) ~ 4.1Td(l — f/2) + 5.1T c //2. 

Including the effect of scatterings we also solved equations (IT41) and ([15]) . and evaluated 
Td and T c . Figure [5] shows the temperatures in the cases with M = lM s _1 . The effect of 
neutrino electron scattering is negligible and the behavior of temperature is well described 
by equations (118!) and ( |T9i) . 

The emerged neutrino spectra with neutrino electron scatterings and pair process are 
shown in Figure [6] (left). Note that the effect of scatterings is so small that scattered 
neutrinos make a small hump at the electron energy, which is overlaid by the spectrum 
of pair process. In the case where / ~ 1 half of the energy dissipated in the corona is 
reprocessed in the disk as thermal neutrinos with lower temperature and the half emitted 
by pair process with higher temperature pass through the neutrino-thin corona. Figure [6] 
(right) shows the amplification of neutrino mean energy, including the scattering and pair 
process. If we take H/r = 0.1 — 1, the mean energy of neutrinos is enhanced by a factor of 
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about 1.5. A thinner corona becomes hotter and emits neutrinos with higher energy because 
the neutrino emission occurs in a small dissipation region as the cooling process. 

4.3. M Dependence 

There is a hope that T c may increase with an increase of M so that the energy amplifi- 
cation factor could increase as M increases. Therefore, we also calculated the amplification 
of the mean energy of neutrinos in the cases where M = O.lM s _1 and 10M Q s _1 . From 
Figure [7] we see that the mean energy of neutrinos from the corona with high accretion rate 
tends to be slightly enhanced since the M dependence of the coronal temperature is weak. 
In the case with M = O.lM s _1 the mean energy of neutrinos is enhanced by a factor two 
with H = O.lr but the neutrino luminosity (and also the absolute value of mean energy) is 
small. On the other hand when M = 10M Q s _1 , the mean energy is enhanced only by 1.2 
even though we neglect the effect of neutrino trapping. If we consider the neutrino trapping, 
the enhancement will be smaller. Thus, there is practically no improvement in the neutrino 
energy, even if M is larger. 

5. SUMMARY AND DISCUSSION 

5.1. Brief Summary 

Magnetically heated corona above the accretion disk may emit non-thermal, high energy 
neutrinos, which, in principle, would lead to an enhancement of the energy deposition rate 
by neutrino annihilation. The improvement is, however, only by a factor of two or so when 
H/r > 1CT 3 . Hence, the energy deposition might still be insufficient to energize GRBs unless 
H/r is extremely small. We see that higher energy of neutrinos is expected if the corona 
form in a thiner region compared to the disk scale height. 

We also see that the effect of the scattering between electrons and neutrinos is negligible 
compared with the cooling by pair process. The ratio of the energy loss rate of the corona 
by electron-neutrino scatterings to the emission rate from the corona by pair process is 
approximately 

(7aT d 4 /8) x (T c /T d )r ^ ^ (T, 

Then, the effect of scattering is negligible at the higher coronal temperature, which is the 
situation which we are considering. We can easily understand this ratio by considering the 
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temperature dependences of the two neutrino processes. The number density of coronal 
electrons is roughly n e oc T c 3 and the cross section for pair process is a e - e +^ uP oc T c 2 , thus 
the emissivity of neutrino pair production is proportional to T c a e - e +^ p n c 2 oc T c 9 . On the 
other hand, the number density of thermal neutrinos is n u oc T^ 3 and the cross section for 
electron neutrino scattering is a cv oc T^T C . Hence, the emissivity of neutrinos scattered by 
coronal electrons is proportional to T c a eu n e n u oc Td 4 T c 5 , where we assume that the energy 
of scattered neutrino is roughly T c . From these simple estimations, we understand that 
the emission of pair process dominates that of electron neutrino scatterings by a factor of 

(r c /T d y. 

From Figure we see that the temperature of the disk does not decrease so much even 
if / increases. This is can be understood in the following way: since about a half of energy 
of emitted neutrinos in the corona is absorbed and heats the disk, more than a half of the 
gravitational energy, Q + , is emitted as thermal neutrinos from the disk. Hence the disk 
temperature decreases by only a factor of 0.5 1 / 4 ~ 0.84 even if / ~ 1 (the case in which most 
of the gravitational energy is dissipated in the corona), compared with the the cases without 
disk coronae. 

We also investigated the dependence of the mean energy of neutrinos on the accretion 
rates. The disk temperature Td,o is proportional to M 1 / 4 since the dissipative energy is 
proportional to M, while the coronal temperature is proportional to M 1//9 if the pair pro- 
cess dominates. Hence the ratio T c /T&q is proportional to M _5//36 . Though the neutrino 
luminosity is proportional to M, even if we neglect the advective energy by neutrino trap- 
ping in the disk, the enhancement of mean energy of neutrinos is small at high accretion 
rates. Thus we cannot expect significant enhancement of neutrino heating by high energy 
neutrinos even at high accretion rates. If we take into account of neutrino trapping, we need 
replace the fraction of energy dissipated in the disk, 1 — /, by 1 — / — / a d v where / a d v is 
the fraction of advective energy. We could neglect the advective energy when / ~ 1 since 
the energy is transported via magnetic process to the surface of the disk on the dynamical 
timescale, which is shorter than the neutrino diffusion time (Ohsuga et al 2002). This will 
lead to improvement of the result of DPN even though the enhancement of the mean energy 
of neutrinos may still be small. 

5.2. Confinement of Coronal Plasma 

The corona might be formed by the energy release of magnetic fields in a thin region 
above the disk. What can then confine the magnetic corona? If we consider the hydrostatic 
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balance in the corona, the thickness of the corona becomes 



V 3 \ r. 

where we assumed the sound speed in the corona as c s = Hence the corona cannot 

be gravitationally bound but expands, unless the radial position of the corona is very close 
to the Schwarzschild radius. The ratio of the cooling time, t coo i ~ llaT c 4 /4g e - e +^ J/P , to the 
free expansion time of the corona, t exp ~ H/c s , is 

*cooi -22 _ 1> ^ (22) 



^exp 



VI' 



where r vV = g e - e +^pif/4(7/8)<7T c 4 is the neutrino optical depth of the corona for pair 
process (RS). Thus, the corona should expand before the coronal plasma is cooled by neutrino 
emission since we now consider the neutrino-thin hot corona, and some processes are needed 
to confine the coronal plasma in a thin region. 

The magnetic pressure could be amplified up to the equipartition value of the gas pres- 
sure in the disk and the amplified magnetic fields will be lifted up in the corona (RS). Such 
strong magnetic fields in the corona might be able to confine the coronal plasma if the 
magnetic pressure is higher than that of relativistic particles in the corona. Since the gas 
pressure in the disk is roughly P gas ~ 10 30 ergs cm -3 for M = lM s _1 (DPN), the magnetic 
fields, whose pressure is comparable to P gas , could confine the corona if T c < lOMeV, i.e. 
H > lCT 3 r in our model (see Fig. [5]). Our estimates above include some uncertainties, e.g. 
the detailed structure of magnetic fields, and we need further studies to see whether such 
magnetic fields really can confine the coronal plasma. 



5.3. Baryon Contamination in the Corona 

Finally, we discuss on the validity of the assumption that the corona consists of pure 
pair plasma with neglegible baryons. In the existence of baryons, neutrino emission by 
the pair capture process may also cool the corona. The cooling rate by the pair capture 
process dominates over that by the pair process if the corona contains sufficient amount 
of baryons, p w > 0.6T 3 10 , where p w is the density of baryons in the unit of 10 10 g cm -3 
and T c l0 = T c /10MeV (Popham, Woosley, & Fryer 1999). Thus, for the validity of the 
assumptions in our model, the baryon density in the corona should be smaller than that in 
the dense disk, whose baryon density is typically 10 10 ~ 11 g cm -3 for neutrino thick disk (e.g. 
Kawanaka & Mineshige 2007). Since the coronal temperature is T c l0 < 1 with a moderate 
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thickness of the corona in our model, H > 10~ 3 r, the baryon density in the corona should 
be less than that in the disk by one or two orders of magnitude. This condition could be 
marginally justified. 

Moreover the distribution function of electrons in the corona could be deformed from 
the Fermi-Dirac type of zero chemical potential. In such a case the neutrino emission by the 
pair process could be suppressed (Kohri & Mineshige 2002). Since the chemical potential of 
electrons is determined by the amount of baryons in the corona, this calculation is beyond 
the scope of this work and should be done in future works. 

The authors are grateful to S. Nagataki and Y. Masada for discussions and comments. 
The numerical calculations were carried out at YITP in Kyoto University. This work is 
supported by a Grant-in-Aid for the 21st Century COE 'Centre for Diversity and Universality 
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Energy Deposition by Neutrinos 
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Fig. 1. — The corona model in hypercritical accretion flows. 
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Fig. 2. — Energy balance in the corona and the disk by neutrino processes. A fraction of 
liberated energy, fQ + , is dissipated in the corona, and (1 — f)Q + in the disk. 
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Fig. 3. — Temperatures of the corona, T c (f tiled), and the disk, T^(open), as functions of /. 
The parameters are M = lM s _1 and the various coronal thicknesses, H/r—1 (circle), 10 _1 
(square), 1CT 2 (upward triangle), 10~ 3 (downward triangle) and 1CT 6 (diamond). Note 
that T c and are normalized by T^o = 4.2 MeV. 
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Fig. 4. — Left: Neutrino spectra emerged from the corona at r = 4r s . Spectra of 
thermal neutrinos with T^o (dot line) and coronal neutrino spectra with (H/r, f) = (l, 
0.91) (short dashed line), (H/r, f) = (l0~ 3 , 0.97) (dashed line) and (H/r, /) = (10~ 6 , 0.84) 
(solid line), respectively. Right: Same as Figj3]but for the mean energy of emerged neutrinos 
normalized by (E)q ~ 17 MeV, the value of the case with / = (no corona). 
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Fig. 5. — Same as Figj3]but for the cases where the neutrino pair process is also considered. 
The solid line shows the temperatures where the neutrino-electron scatterings are neglected. 
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Fig. 6. — Left: Same as FigJH but for the case with neutrino pair production. The param- 
eters are (H/r, /) = (1, 0.96) (short dashed line), (H/r, f) = (10~ 3 , 0.90) (dashed line) and 
(H/r, f) = (10~ 6 , 0.94) (solid line). Right: The amplification factor of the neutrino mean 
energy, same as Figj3] 
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Fig. 7. — Same as Figj6] (right) but for the case with M = O.lM s _1 (left) and M = 
lOM s _1 (right). In the case with M = 0.1M Q s _1 the neutrino mean energy is normalized 
by 9.7 MeV, and with M = lOM s _1 normalized energy is 31 MeV. 



